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Abstract 

We compute the renormalized trajectory of ^^-theory by perturbation theory 
in a running coupling. We introduce an iterative scheme without reference to a 
bare action. The expansion is proved to be finite to every order of perturbation 
theory. 



1 Introduction 



In Wilson's renormalization group |[W71| , [WK74|| ultraviolet and infrared limit stand 
for an infinite iteration of block spin transformations. Consider for instance the ultravi- 
olet limit of an asymptotically free model at weak coupling. There the point is to keep 
couplings under control which grow under a block spin transformation. Such couplings 
are called relevant. In weakly coupled models they can be identified by power counting. 
Renormalization of a bare action amounts to sending it through an increasing number 
of block spin transformations. The image is the renormalized action. For this limit 
to exist the bare couplings need to be tuned as the number of block spin transfor- 
mations is increased. This renormalization scheme has been beautifully implemented 
both within and beyond perturbation theory. Let us mention the work of Gallavotti 
p85 |, Gawedzki and Kupiainen [pK77 , |GK84 |, and Polchinski | P84 | as a guide to the 
extensive literature. The renormalized actions are located on a low dimensional curve, 
parametrized by renormalized couplings. In the case of 0^-theory we can consider a 
one-dimensional subspace and speak of a renormalized trajectory. The idea is to ar- 
range the renormalization group flow generated from the bare action as to converge 
to this renormalized trajectory in the process of renormalization. The renormalized 
trajectory of an asymptotically free model is pictured as unstable manifold of a trivial 
fixed point. 

Although this picture has been behind block spin renormalization since the very 
beginning it has not yet been formalized to an approach free of bare action. This 
paper is a contribution to fill this gap. It extends the analysis begun in |[WX94| , |RW95| 



for the hierarchical approximation to the 0^-model with discrete momentum space 
renormalization group. The renormalized trajectory is here defined as a curve in the 
space of effective actions which passes through the trivial fixed point and whose tangent 
at the trivial fixed point is a ^'^-vertex. The dynamical principle which proves to 
be sufficiently strong to determine this curve at least to all orders of perturbation 
theory is stability under the renormalization group. With stability we mean here 
that the curve is left invariant as a set in the space of effective actions under a block 
spin transformation. A renormalized action always comes together with a sequence 
of descendants generated by further block spin transformations. Even in the case of a 
discrete renormalization group this sequence proves to consist of points on a continuous 
curve in the space of effective actions which is stable under a block spin transformation. 
It is the computation of this curve in a vicinity of the trivial fixed point we address. 

Given a block spin transformation, we may distinguish between the following dif- 
ferent renormalization problems. The first problem is an initial value problem, the 
analysis of the renormalization group flow started from a particular bare action. The 
second problem is a mixed boundary value problem, where the relevant parameters are 
prescribed on a lower scale, the irrelevant parameters on a higher one. The first prob- 
lem is appropriate for the infrared limit of a Euclidean field. The second is appropriate 
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for its ultraviolet limit. The third problem is the question of fixed points. In this paper 
we consider a generalization of the fixed point problem. We will look for interactions 
which remain invariant up to a (one dimensional) flow of a coupling parameter. A 
requirement of finiteness will substitute for boundary data. Of course the problems are 
interrelated. In particular the mixed boundary value problem is a method to obtain a 
solution of our generalized fixed point problem in a scaling limit. 

The result is an iterative form of perturbation theory in a running coupling. Its 
closest relative in the literature is the renormalized tree expansion of Gallavotti and 



collaborators ||G85| , pN85|| . Our expansion will however not be organized in trees. 
Furthermore, it will from the very beginnining be free of divergencies piled up in 
standard perturbation theory. Surprisingly it will allow to treat relevant and irrelevant 
couplings on the same footing. It will involve neither bare couplings nor renormalization 
conditions in the original sense. The expansion will be presented for the 0^-trajectory 
in four dimensions. Most of the analysis, in fact everything except for the treatment 
of the wave function term, works in arbitrary dimensions. We therfore leave D as a 
parameter in the equations. Thereby dimension dependence of scale factors is exhibited. 
The three dimensional case requires a modification which will be explained elsewhere. 
See ||RW95|] for a treatment of its hierarchical approximation.^ 

To make contact with the physical world one has to supply one more piece of 
information. One has to assign a scale to a point on the renormalized trajectory. In 
the presentation of this paper we will maintain a unit scale throughout the computation. 
In view of asymptotic freedom the four dimensional trajectory is suited for the infrared 
limit of massless 0^-theory at positive coupling. It is also suited for the ultraviolet 
limit of massless 0^-theory at negative coupling, as been promoted by Gawedzki and 
Kupiainen ||GK85|| . In this paper we will restrict our attention to the effective action. 
Its relation with Schwinger functions is discussed for instance in PG95|| . The iterative 
solution of our equations can also be viewed as an improvement program organized in 
powers of a running coupling. 



2 Renormalization group 

The below analysis will be done in terms of a discrete momentum space renormalization 
group transformation. A number of applications of which is discussed for instance 
in recent lectures PG95|| by Benfatto and Gallavotti. Let us consider the following 
discrete block spin transformation TZ on some space of interactions V{(j)) of a real 
scalar field (j) on Euclidean space IR^. Let TZ be composed of a Gaussian fluctuation 
integral with covariance F and mean ip with a dilatation S of the background field if). 

^In three dimensions the expansion parameter has a non zero scaUng dimension. As a consequence 
the flow of the mass term shows up a second order correction proportional to the logarithm of the 
running coupling. This type of corrections requires of a double expansion in the running coupling and 
its logarithm. 



2 



Let the fluctuation covariance be defined by 



Hp) = -2 - x{Lp)) , (1) 

where x{p) is a momentum space cutoff function. Its purpose is to make T{p) decrease 
fast outside a momentum shce L^^ < \p\ < 1 set by a scale parameter L > 1. A 
convenient choice is the exponential cutoff 

X(p)=e-^'. (2) 

It will be used in the following. Other choices however work as well, for instance Pauli- 
Villars regularization. Then (|I]) defines a positive operator on the subspace of L2(1R'^) 
consisting of functions /(x) with zero mode /(O) = 0. Let d/ir(C) be the associated 
Gaussian measure on field space. Recall its basic property 

d/ir(C)e(^'^) = e^(^'^^) (3) 

and consult Glimm and Jaffe [pJ87|| for further information. Let the fluctuation inte- 



gral then be given by the average of the Boltzmann factor Z{(j)) = exp{—V{(j))) with 
respect to d/ir(C), shifted by an external background field ip- Let us introduce the 
notation 

{Z)^^^ = JdMC)zi^ + C) (4) 

for this average. We can think of the momentum slice L^^ < |p| < 1 as a portion of 
momentum space degrees of freedom which is integrated out. The integration of another 
portion is prepared for by a dilatation of the background field. Let this dilatation be 
given by 

S^ix) = L'-^^P (I) . (5) 

The exponent a = 1 — D/2 is the scaling dimension of a free massless scalar field. 
Anomalous rescaling will not be considered here. Non-anomalous rescaling applies 
(at least) to small perturbations of a free massless field. The renormalization group 
transformation is then defined by composed with (j^). The following analysis will 
be done in terms of the potential V{(j)). The method will be perturbation theory. The 
matter of stability bounds on Z{(j)) will not be addressed here. The renormalization 
group transformation for the potential then reads 

W(^) = -log((exp(-l^))r^J. (6) 

We will restrict our attention to even potentials V {—</)) = V{(f)). The transformation 
preserves this property. Potentials differing by a field independent constant will 
be identified. V{(j)) can for instance be normalized such that V^(0) = 0. To maintain 
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normalization, (^) then should be supplemented by a subtraction of TZV{0). Technically 
this constant is proportional to the volume, infinite in infinite volume. Eq. (j^) therefore 
requires an intermediate volume cutoff to make sense. We will wipe this technicality 
under the carpet, keep (H) as it stands, and ignore field independent constants. When 
we perform a renormalization group transformation we will speak of V{4>) as bare and 
of 7lV{(j)) as effective or renormalized potential. It should however be kept in mind 
that only degrees of freedom in one momentum slice are integrated out in a single 
renormalization group step. ( will be called fiuctuation field and ip background or 
block spin field. The term potential is sometimes reserved for local interactions. Here 
potential will be used synonymous with full interaction including nonlocal interactions 
generated by the renormalization group. The block scale L will be kept fixed in the 
following. It should not be confused with a full momentum space cutoff of a Euclidean 
field. A typical value of L is two. A full cutoff could be an iV'th power of L. 



3 Trivial fixed point 

The renormalization group transformation has a trivial fixed point K(0) = 0. This 
fixed point is the free massless scalar field. Eq. has in fact been designed upon a 
momentum space decomposition of a free massless field. The linearized renormalization 
group transformation at this fixed point is given by a Gaussian expectation value 

Bv^nOi^) = 5^ , (7) 

shifted by a rescaled external field. It is diagonalizable. The eigenvectors are normal 
ordered products. We will represent the potential in terms of normal ordered products. 
Let us therefore briefiy recall some basic facts about normal ordering. Normal ordered 
products with normal ordering covariance v are generated by 

The linearized renormalization group for this generating function is an exercise in 
Gaussian integration (|]). The result is 

/ d/ir(C) : e(^^+^'^) :„=: e^^'^^^) :r. • (9) 
The generating function (H) is preserved up to a (transposed) dilatation 

S^f{x) = L'^^l'f{Lx) (10) 
of the source and a linear transformation 

rv = S-\v-T){S-^f (11) 
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of the normal ordering covariance. This hnear transformation generates a flow of 
normal ordering. It can be thought of as a residual renormalization group flow taking 
place besides more interesting dynamical effects. It has a line of fixed points 

Hp) = \ ixip) - c) (12) 

parametrized by C. Let us select the point C = as normal ordering covariance in the 
following, a massless covariance with unit ultraviolet cutoff. Since it remains invariant 
under ([TTD it can be safely suppressed in the notation. It follows that 

= j d^xi ■ ■ ■ d^x„ : ^(xi)'"^ ■ ■ ■ 0(x„)"^" : (13) 

is an eigenvector of the linearized renormalization group (|^ with eigenvalue L'^ . The 



exponent is a = nD + (mi + • • • + m„)(l — -D/2), the scaling dimension of (0). A 
prominent member of this family is the 0^-vertex 

(9(0) = y d^x : (P{xf : (14) 

with scaling dimension A — D. It is therefore called relevant in D < 4, marginal 
in D = 4, and irrelevant in D > 4 dimensions. General eigenvectors are given by 
homogeneous kernels in real or momentum space and also involve derivatives of fields. 
See the review |W76| ] by Wegener. 



4 Perturbation theory 



In a vicinity of the trivial fixed point the transformation (^ can be computed by means 
of perturbation theory. The perturbation expansion for the effective potential reads 



^1 



T 



(15) 



The superscript ^ indicates truncated expectation values. Truncated expectation values 
are defined by 



-1 , 



\i=l 



(16) 



J Ai = 



=A„=0 



Notice that the truncated expectation values ([^;]")r5^ the cumulants associated 
with the moments (^'^)r5^- The perturbation expansion (|T^) is therefore also known 
as cumulant expansion. Let us consider the situation when the potential comes in form 
of a power series 

oo n 



n=l 



(17) 
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of a coupling parameter g. The zeroth order interaction will always be V^^\(f)) = in 
the following. In the bare perturbation expansion for a single renormalization group 
transformation the effective interaction is expanded again in the bare coupling. Let us 
introduce the notation 



oo n 

nv{ij\g) = Y.^-TmY''\^) (18) 

for this expansion. The individual orders of perturbation theory are given by sums of 
truncated expectation values 

invrw^t i^ j: )(^^'^^;- .v^^-^):,,. 

/i,...,/^elN ^ ^ 

h-\ Vlm = n 

(19) 

The sum over (Zi, . . . , 1^) is restricted to m-tupels in {1, ... , n}™ such that + - ■ ■+lm = 
n and is finite. The multinomial coefficient is given by 

(20) 












(KV)W(^) = 








(KV)<"(« = 






3/y(2) y(i)\^ +/V^(i) 



with = n — li — ■ ■ ■ Im-i- To lowest orders we find the following explicit expressions 
for the effective interactions 

(21) 
(22) 
(23) 

Let us have a closer look at two contributions to (|1^). The highest order bare interaction 
appearing to order n is V^^\(j)). It contributes a term (v'^^'''^^ It is thus transformed 

according to the linearized renormalization group (0). The perturbative corrections to 
the linearized renormalization group depend on lower orders V^"^\(f)) with 1 < m < 

n — 1 only. The first order V^^\(j)) contributes by itself a term (—1)"^^ ([^*^^^' ]")r 5i/) 
to (JZV)^"'\'il') . These interactions need to be carried along immediately to order n 
when some first order interaction enters the game. In a minimal scheme no further 
interactions would be introduced to this order. 

(|T9|) leaves us with an expansion for the effective potential in terms of the bare 
coupling. This expansion is not appropriate for an iteration of renormalization group 
transformations. The appropriate expansion is an expansion in powers of the effective 
coupling. This requires a reorganization of (|l^) in powers of the effective coupling. Let 
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this effective coupling come as a power series 



oo n 

m = E '-jbn (24) 

n=l "'• 

in the bare couphng. The function f]{g) is not a Callan-Symanzik /3-function in a hteral 
sense but its block spin analogue. The inverse reorganization goes as follows. Let the 
effective potential be given as a power series (0) of the effective coupling. Then we 
can substitute (|2^) and expand the result in powers 

oo n 
n=l 



of the bare coupling. The coefficients in this reorganized expansion (|25|) are given by 
{Vo(3)(-){ij)=±l^ E ( , ""j ]bH---KV'-"^\^). (26) 



m=l 



G IN 

h + \-lm = n 



in terms of the coefficients of ( p^TD and (0). The zeroth order of (p^ ) is put to bo = 0. 
To lowest orders (|2^) is given by 



{Vopy'\^) = 6i\/«(V^), (27) 
(Vo/5)(2)(^) = 62V^«(V;) + (6i)V(2)(V;), (28) 
(\/o/5)(3)(^) = b3V^'\ij) + 3hb2V^^\^) + {bifV^^\ij). (29) 

This combinatorial exercise completes the setup for perturbation theory. 

If we would iterate (0) we would end up with the tree expansion of Gallavotti 
|p85|| . See ||1:''HRW88| for a pedagogical account of this well developed technology. 



We will not organize perturbation theory in terms of trees. Also we will not start 
with a bare perturbation expansion but directly attack the renormalized series. Let us 
conclude this section with the remark that the effective potential (plSj) is the generating 
function of free propagator amputated connected Green's functions with vertices ( [T7| ) 
and propagator F. The terms in ([19D can therefore be expressed in terms of connected 
Feynman diagrams. 



5 0^-trajectory 

Let us pose the following renormalization problem. We seek a potential V{(j)\g) de- 
pending on a running coupling g with the following properties, ("running" will be 
explained below.) 
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1) V{(j)\g) is a power series 

CO n 

n« = E^^^"n0) (30) 



n=l 



in the coupling parameter g. 

The zeroth order is here V^^^cj)) = 0. We will treat of ( pOD as a formal perturbation 
of the trivial fixed point. The important question of summability of (|30D will not be 
addressed. 

2) The first order in (p^j is a (p'^-vertex 

V^'\^) = lJd''x:<P{xy:. (31) 



The n 'th order in ( [3^ j is a polynomial 



n+l 1 „ 

V^^\<P) = ^ / dxi . . . dX2mV2t\xu . . . , X2™) : 0(xi) ■ ■ ■ <j){x2m) : (32) 

iii tiie field 0. 

Let us imagine a bare theory consisting purely of a first order vertex (|3TD, possibly 
shouldered by second order mass and wave function counterterms. The effective in- 
teractions generated in the course of its renormalization group flow will not remain 
of this simple form. The first renormalization group step already leaves us with an 
infinite set of higher order vertices, among which are for example second order mass 
and wave function terms. A main theme of this paper is to add all these higher order 
interactions to the bare theory from the beginning since they will anyway be generated 
once (|3TD enters the theory. In a minimal scheme no other vertices are introduced than 



those enforced by the presence of (0). We can then think of (0) as the germ of the 
theory. Other trajectories emerging from the trivial fixed point are of interest as well, 
for instance the 0^-trajectory. There the first order is a 0^-vertex. Our leitmotiv here is 
to keep the effective interaction as minimal as possible. The highest connected vertex 
generated from n first order 0^- vertices has 2n + 2 legs. Higher vertices will appear at 
this order only if they are introduced by hand. excludes this option. This form 
of potential iterates through the renormalization group. Field independent terms are 
discarded. We consider even powers of fields only. 

4) The kernels in ( P^ ) are Euclidean invariant distributions. They are given by Fourier 
integrals 



V^'Jixi, . . .,X2m) = J -j—^ ■ ■ ■ -^^^e ^'=1^' \2tx) 5\^pij ^^(pi, . . . ,P2. 



(33) 
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With the 6-function removed their Fourier transforms are symmetric C°° -functions on 
momentum space H 



D 



X 



X JR^ . They satisfy the bounds 



P 



dp" 



f>(") 

^2m 



sup 

2n 



P"^^m(Pl,---,P2m) 



-e(|pi| + ...+ |p2m|) 



< oo. 



for any e > 0. Here A4 2m denotes the subset of momenta {p 



1) 



^P2r. 



G IR^X' 



(34) 



■xH 



D 



with pi + 



+ P2m = 0. Furthermore, 



a 



N^™-^ is a multi-index and 



The kernels in (|33|) are the unknowns in this approach. (0) is meant as requirement on 
what kind of kernels we will hold look out for. This condition of finiteness substitutes 
for boundary data which is necessary for instance in |p85 | and | P84 |. On the practical 
side we want to be certain that our perturbation theory is finite to every order. The 
Loo^e-norm in momentum space is a convenient but not the only possible criterion. 
See ||FHRW88| for further inspiration. The large momentum bound implied by ( p4|) is 
rather wasteful. The accurate large momentum behavior is polynomial in powers and 
logarithms of momenta. Their origin is a Taylor expansion of the quadratic and quartic 
kernel in four dimensions. The Taylor expansion is compelled by power counting. The 
remainders are irrelevant. The norm estimates present bounds on these irrelevant 
remainder terms. The Taylor coefficients are of course also required to be finite. This 
is understood when we speak of smooth functions on momentum space. Smoothness is 
a little more than what is needed. Existence of required derivatives at zero momentum 
(or some general subtraction point) together with ,;-bounds on the remainders is 
fully sufficient. In the ultraviolet problem with exponential cutoff we can afford the 
luxury of smoothness. Euclidean invariance reduces the Taylor coefficients in four 
dimensions to a mass term, a wave function term, and the 0^-vertex. Note that @ 
respects Euclidean invariance. It is therefore broken only if we break it by hand. Note 
also that translation invariance by itself affects momentum space power counting. 



5) The four point kernel ( pSj ) is 



vt\o,o,o,o) = Sn, 



(35) 



at zero momentum. 
The zero momentum condition 



I) is part of the definition of the expansion parameter 
g. If g is traded for another parameter g'{g) = g + 0{g'^) then ( p5D is not true anymore. 
There is nothing wrong with a redefinition of g, although it looks unnecessary to 
introduce a 0^-vertex at some higher order when it is already present to first order. 



6) There exists a power series 



oo n 
n=l "-• 



(36) 
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such that 

nViij\g) = Vii;\Pig)) (37) 

to every order in g. 

This condition is the core of our approach. It says that V{(l)\g) remains of the same 
form under the renormahzation group up to flow of g. It is therefore called running 
coupling. ( p7D will prove to be strong enough to determine both V{(f)\g) and P{g) to 
every order in ^f. In D = 4 dimensions this scheme will require a modification due to 
wave function renormahzation. We will need to introduce a first order wave function 
term. But let us postpone this modification for the moment. A potential with the 
property (^) will be said to scale. 

It is amusing to think of the solution V^(0|-) : g ^ y{4>\g) as a parametrized curve in 
interaction space. Anticipating a future renormahzation group geometry we can then 
say that: 1) V^(0|-) visits the trivial fixed point a.t g = 0. 2) The tangent to ^(0|-) at 
g = is the 0^-vertex. 3) ^((/>|-) is stable under the renormahzation group as a set in 
interaction space. We call this curve 0^-trajectory. 



6 Scaling equations 



Expanding both sides of (|37|) in powers oi g, we obtain a system of scaling equations 
for the perturbations series (|30| ) and (|36D. This system is organized into a recursion 
relation. From (pT]), (p2D , (|23|) and (|27|) , (|28|) , (^) we find to lowest orders equations 
of the form 

(\/«)^_^^- 6i\/«(^) = 0, (38) 
{v^\^^-{hrv^'\^P) = 6,V«(^) + (V«;V«)^^^^, (39) 

(40) 

To first order, scaling requires V^^\ip) to be an eigenvector of the linearized renormal- 
ization group, and bi to be the eigenvalue. This is the case for the 0^-vertex with the 
familiar eigenvalue bi = L^~^ . The first order equation is special in that it is homoge- 
neous. To higher orders we meet inhomogeneous linear equations. Their general form 
for n > 2 is 

V-W)^ (6,)"rW(^) = 



J ^ n 



/i,...,/^GlN ^ 

h-\ Vlra = ri 



10 



n— 1 1 y T \ 

T.^^\bh ■ ■■buV<^'-\^P) + (-l)'" (^(^^); ■ ■ ■ ; J- (41) 

We give the name L^^'^~^^ K^'^\'ijj) to the right hand side of (|4l|) . It is determined by 
the effective potential V^'^\ip) to lower orders 1 < m < n — 1. Eq. (|1) poses an 
inhomogeneous linear problem for V^"'^ (ip) . Its left hand side is diagonalized by normal 
ordering. Insert the expansion (^) in normal ordered powers of fields to obtain 

n+l „ 

Lrr.(D+2)y(n)^^^^^ . • . , Lx^m) ' L^^^"^) V ^ (xi , . . . , X,j) . (42) 



The linearized renormalization group thus performs a scale transformation of a kernel.^ 
The difference of exponents defines an order dependent real space power counting 

a{m,n) = m{D + 2)-n{A-D). (43) 

The right hand side of (|4l]) can now be expanded in the same manner. We introduce 
the notation 

n+l 1 „ 

irW(^) = J2 TTT^ / d''xi---d''x2m:Hxi)---^{x2m):Ki;^{xi,...,X2m) (44) 

for this expansion in normal ordered powers of fields. From (^) and ( ^^ together we 
then deduce the equations 

L-My^n)^^^^^ . . . , LX2^) - VS{Xl, • • • , X2m) = K^£.{Xl, • • • , X^m) ■ (45) 

The problem of renormalized perturbation theory has thus been reduced to the solution 
of (1^). Recall that the right hand side is determined by vertices of lower orders and 
cutoff propagators. The kernels are distributions on copies of real space. Fourier 
transformation turns (|45|) into 

^^Kn)^W(|^ . . . , ^) - VSiPl, . . . .V2ra) = K^^{Pl, • • • ,^2^)- (46) 

Here the 5-function due to translation invariance has been removed. Eq. (|46|) is the 
main dynamical equation in this approach. The task is to search for finite solutions in 
the sense of (|3^). In momentum space we find an order dependent power counting 

a{m,n) = D -m{D -2) -n{A- D). (47) 



^Homogeneous kernels therefore give eigenvectors of the linearized transformation. 
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In the sequel the attributes relevant, marginal, and irrelevant will be used to distinguish 
whether (^) is positive, zero, or negative. Notice that ( ^Tj) becomes 



m 



D 



4 - 2m, D = 4, 



(48) 
(49) 



in three and four dimensions respectively. Thus in three dimensions three interac- 
tionsQ are non- irrelevant, while there are infinitely many in four dimensions. Labelling 
interactions by a pair (m, n) we have: 





(m,n) 


relevant 


marginal 


D 








3 


(1,1) 


(1,2), (2,1) 


4 


(l,n) 


(2,n) 



The conclusion is of course that the theory is super-renormalizable in three and renor- 
malizable in four dimensions. In the BPHZ scheme this leads to the complication that 
the renormalizable case has infinitely many divergent graphs. In the present approach 
(with no divergent graphs) the renormalizable case will not be more complicated than 
the super-renormalizable one. 

Consider a general difference equation of the form 



( f ) - fip) = 9ip) 



(50) 



for functions on IR^. Assume that L > 1, a is integer valued, and g G C°°(1R^). We 
distinguish two cases, the irrelevant case with a < and the relevant case with a > 0. 

Let a < and g E C°°(IR"). 

1 ) The series 

P 



fip) 



(51) 



is uniformly convergent on compact subsets of IR .It deGnes a function f G C°°(1R"'). 

2) The function f given by ( 15ID is a solution to (]5^j. It is unique in the space of smooth 
functions. 

3) Let \g\oo,e < Then f satisfies the bound 

1 



l/lc 



< 



1^1 



(52) 



^Interaction here refers to the fuU kernel and comprises its value as well as derivatives at zero 
momentum. 

"Recall that |5|oo,e = suppgjRW \g{p)\e:cp{-e\p\). 
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If the derivatives of g are Loo,e-bounded, so are the derivatives of f. 

Formula (|5TD is obtained by iteration (0). Its convergence follows from the Weierstrass 
condition for uniform convergence |[WW86|| . It is checked to satisfy (0). The difference 
of two solutions to (50) satisfies the homogeneous equation L^f (f ) - f{p) = 0. When 
a is negative its only smooth solution is zero. The Loo,e-bound is elementary. 

Thus eg . (^OD has been solved in the irrelevant case. In the relevant case the trick 
is to reduce the degree of relevancy by taking derivatives. Deriving (^D|) once, we have 



dp \Lj dp dp 

A partial derivative therefore reduces powercounting by one unit. The value of the 
derivative at the origin is obtained by evaluation of (|53D . The strategy is therefore 
a Taylor expansion with remainder term to an order where the remainder becomes 
irrelevant. 

Let (J > and g e C°°(1R"). Then 

Sip) = E^;^!^^^^^ E^_^i'*(i-<r|^Mp" (54) 



is a solution to ( l5^ ) if and only if 



(i"-'-' - 1) ^(0) = -g^iO). \a\<a. (55) 

i.-H*^m-^(rt . ^(rt, + (56) 

dp"' \LJ dp'^ dp"' ^ ' ' ' ^ ' 



The proof is obvious. As a consequence smooth solutions exist only when all marginal 
Taylor coefficients of g vanish at the origin. The marginal Taylor coefficients of the 
solutions are free parameters. 

Let (T > and g E C°°(1R"). 

1 ) Smooth solutions to ( ^ ) exist if and only if 



dx'^ 

for all multi-indices with la I = a. 



-(0) = (57) 



2) In this case, we obtain a set of solutions which can be parametrized by the Taylor 
coefficients 

c. = ^(0) (58) 
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for all multi- indices with \a\ = a. 

The proof is obvious. The grain of sah is that the marginal Taylor coefficients are 
not determined by eq. (|50D . The relevant Taylor coefficients and the irrelevant Taylor 
remainder are however explicitly known from (|5TD, (pSj), and (^). 

If the Taylor expansion (0) is pushed farther than to order a, then also the Tay- 
lor coefficients with negative power counting obey (|55|). The iterative and the direct 



solution of (55) coincides when a — \a\ < 0. Thus in terms of Taylor coefficients there 
is surprisingly no difference between the relevant and the irrelevant ones. Finally, re- 
mark that also the situation with non-integer scaling dimension a has applications in 
the renormalization group. It is possible to cook up models with irrational scaling 
dimensions. Then the scheme is particularly simple because marginal coefficients are 
absent. 

With this interlude we are essentially finished. The irrelevant kernels in (|i6| ) are 
solved by iteration. The non-irrelevant kernels are Taylor expanded. The Taylor coef- 
ficients are solved directly. 



7 Marginalia 

Marginal eigenvectors require a separate treatment. The following scheme is designed 
for 0^-theory in D = 4 dimensions. More generally it applies to models whose coupling 
is dimensionless. The underlying idea is that marginal eigenvectors possess a logarith- 
mic flow. If also the coupling is marginal we may express the former logarithmic flow 
in terms of the latter. 0^-theory in D = 3 dimensions for instance calls for a different 
approach which will be explained elsewhere. There the coupling flows powerlike and 
one needs to introduce terms proportional to logarithms of the coupling to handle the 
marginal eigenvectors. For the rest of this section we put D = 4. Let us flrst do the 
computation of the /3-function without wave function renormalization to explain the 
scheme in a slightly simpler setting. It is convenient to write it in terms of projec- 
tors on the marginal eigenvectors. Their meaning is nothing but Taylor expansion in 
momentum space. For a general Euclidean invariant functional 

^(^) = Yl 7r~r / '^^^1 • • • (^^^2mV2miXi, X2m) ■ (Pi^i) ■ ■ ■ ^(Xam) : (59) 

deflne projectors 

P2,2V^(0) = IJ d^xi : (90(xi))2 : ^ J d^X2(xi - X2) V2(xi, X2), (61) 
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on the marginal eigenvectors 



OaM = ^ / d^x : :, (62) 

^2,2(0) = lU^x: {d(t){x)f : . (63) 



2 

These eigenvectors correspond to the particular kernels 



V2n{Xi, . . . ,X2n) = Sn,lS{x2 - Xi)6 {X3 - Xi)S {X4 - Xi) , (64) 

dx'i 



V2n{Xi, . . . ,X2n) = 5n,l Tr^'^(a;i - ^2) , (65) 



respectively, symmetrized in their entries. It is clear that they are reproduced by the 
projectors ( |60D and (|6T|). Furthermore it is clear that the projectors satisfy 

P4,o(V^)r,5^ = L^-^'PifiVi^P), (66) 
P2,2(l^)r,5^ = P2,2V^(^). (67) 

The wave function term (^) is marginal (at order zero) in all dimensions and the (f)'^- 
vertex ( |62D in four dimensions. (|6TD and (|60|) are projectors on the marginal content 
of (|59D . Let us now consider the /3-function. For all orders n > 2, the definition of the 
couphng constant implies that 

P4,oV^^"n0) = 0. (68) 

That is the higher orders do not contain an eigenvector ([631) . Applying ([6T| ) to both 
sides of (^) , and using (^) to conclude that the left side projects to zero, we find 



n—l / -| \m+l 

&nP4,0V^«(V^) = E 

e IN 

h + ■ ■ ■ + Im = n 



m=2 



n 

h - ■ ■ Im-l 



P4,o(^^'^V--;^^''"^)J!_5^. 

The projector P4 q selects the contributions built from lower order vertices, fluctuation 
propagators, and normal ordering propagators to the effective n'th order 0^-vertex. 
This equation determines the coefficient 6„. 

Let us now also introduce a wave function renormalization. This requires a slight 
modification of the renormalization scheme used so far. The second order condition on 
the wave function part is 

h2^2,2V^'\^) = -P2,2 (V^(^^ ^^'^>r,5^ ■ (70) 
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From it we see that we cannot leave the first order to remain a pure 0^- vertex ( |3TD unless 
the right side of ( [70D does not contain a wave function term. This is unfortunately not 
the case. A wave function term is indeed generated to second (and arbitrary) order 
of perturbation theory. Therefore we have to include a wave function term already to 
first order in the running coupling and replace (|3ll) by 

V^'\ct^) = \J^''^-- ■■ +1 / : m^)f ■■ ■ (71) 

Here Ci is a new parameter, the first order wave function renormalization. Its value 
is determined by the second order equation (^). On a first inspection this equation 
looks quadratic in Ci- This is not the case since 

P2,2 (02,2; 02,2)ls^ = 0. (72) 

The reason is that the fluctuation propagator (|1|) is regular at zero momentum. Another 
fact which is pleasantly welcome is 

P2,2(O2,2;O4,0)r,5V^ = 0. (73) 

As a consequence the first order wave function term does not alter the second order 
coefficient 62 of the /3-function. A wave function term on an external leg of the 0^- 
vertex makes the resulting amplitude vanish at zero momentum due to the regularity 
of (|1|) at zero momentum. As a consequence Ci does not appear in 

h2P,,oV^'\i^) = -P4,o {V^^: ■ (74) 

from which 62 is computed. Thus we can first find 62 from ([7^) and then (^1 from (|70D . 

Thereafter we can proceed with the computation of the full second order. This 
scheme generalizes to arbitrary order of perturbation theory. Assume that we have 
computed K^™^) (0) for all orders \ <m <n — l except for the wave function renormal- 
ization Cn-i at order n — 1. Then the scaling equation 

hn^,,oV^'\i^) = -nP4,o (^^'^; ^^""'^)r.5^ + (75) 




allows us to determine the n'th order coefficient hn of the /3-function. Analogous to 
the second order case one does not find Cn-i on the right side of ([75|). To make the 
mechanism explicit let us split off the wave function term from the effective interaction 
and write 

V^-^\lP) = Cm02,2m + W^'^\i,). (76) 
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Then by the same reasons as above it follows that 

P..{v'^>-M-'-\^^ = P,,{W^'W^-^>),,^ (77) 

is independent of Cn-i- But all other contributions in ( [75D stem from lower orders only. 
The final piece of gymnastics is to find (n-i from 

n(n — 1] 



2 



n— 1 -j^ 



/i,...,/^ G IN 

h-\ \-lm = n 

.... - nF2,2 {V^'^; V^-''^)l^^ + (78) 

(^(1); . . . ; yd); + (yd); . . . ; yd))^^^^ . 

Again the right side does not depend on the order n — 1 wave function renormalization 
constant (n-i because 

P..(V'";V'"-))^^^^, = P,.(H"');H""-))^^^^. (79) 

Notice that the wave function renormalization constant is hidden notationally on the 



left side of (|7|) in 

P2,2V^''-'\^)=Cn-l02,2{^). (80) 

Notice also that it comes together with pre-factors which do not vanish. The iterative 
scheme is now complete. In summary it goes as follows. To order n we first compute 6„ 
from (^), then Cn-i from ([78|), and thereafter the effective interaction y*^"^(?/') except 
for (n- Then we iterate the computation in the next order. 

8 Estimates 

In the iteration the lower order vertices are convoluted with hard and soft propagators 
to form the next order inhomogeneous scaling terms. In this section we present an 
estimate which shows that the integrals in this process converge. Consequently the 
expansion is finite to every order of perturbation theory. In a nutshell the expansion is 
finite, because the induction step involves perturbation theory with cutoff propagators 
and smooth vertices which do not grow too fast at large momenta. The fluctuation 
propagator has two sided cutoffs. The estimates 

||f||oo,-2. = sup {|f(p)|e2^W}<oo, (81) 



|r||i,-2e = I -^\np)\ e^^H<oo, (82) 
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for e sufficiently small, should therefore cause no surprise. We save a small amount 
of exponential decrease to bound the large momentum growth of the vertices. It is 
then fairly obvious that a perturbation theory with this fluctuation propagator and 
-^oo,e-vertices is finite. The normal ordering propagator has an ultraviolet unit cutoff 
but no infrared cutoff. It satisfies 

_2e < OO. (83) 

But this integrability in dimensions larger than two is sufficient to bound loop integrals 
with normal ordering propagator. Fortunately they never occur at external legs, where 
they would spoil the Taylor coefficients. 

Recall that the induction step to order n involves the computation of 6„, Cn-i; ^iid 
thereafter V^'^\il)) (except for C„). We present the argument for finiteness of V^'^\'il)) to 
some detail. Finiteness of the coefficients then follows by the same reasoning. Let us 
assume bounds on 6„, Cn-i, V^'^~^\tl)), and respective lower orders. The right side of 
the recursion equations (^) consists of two contributions, the dynamical contribution 

" (-Ip ( n (84) 

/i,...,/^GlN ^ ^ 

^1 H Vlm = ri 

from the fluctuation integral, and the kinematical contribution from the flow of the 
couphng. 




m=l 

h + ■ ■ ■ + Im = n 

The second one immediately inherits a bound from the induction hypothesis on the 
lower orders. Recall that it arises from a re-organisation of the interaction as a poly- 
nomial of the bare coupling into a polynomial of the effective coupling. This re- 
organisation is finite because the lower order coefficients of the /3-function are finite. 
To estimate the (P^, we break it down to a sum of Feynman amplitudes. Fortunately 
we do not need more information besides that it can be written as a sum of finitely 
many amplitudes with certain properties. Each of the lower order interactions is a sum 
of vertices V^^^ip) = Eit\ V^i\^) of the form 

vSi^) = J^I d""^! ■ ■ ■d^'xsfc : ^(xi) ■ ■■^P{x2k) : ^?(xi, . . .,X2k). (86) 
We estimate individually each contribution in the truncated expectation value 

k\ — 1 — 1 
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for the sake of notational economy. Insert the expression 
an estimate for 



. Then the task becomes 



T/(^l). . T/('m) 




The truncated expectation value in (|88D contains a product of clusters, each cluster 
being a normal ordered products of fields. The index j defines a colouring of the 
clusters. Eq. (|88| ) is computed in two steps. Step one is 



n 



1=1 



/iC{l,...,2fei} ImC{l,...,2k^} j=l iaij 



X 



n 



iG{l, 



n ax,, 

...,2k,}\I, 



(89) 



The truncated expectation value is here a standard one with mean zero. We sum over 



all subsets Jj of {1, . . . , 2kj}. The interpretation is that for each field 



(X 



belonging 



to cluster j we decide whether it be a rescaled background field Sip{xj^i) or a fluctuation 
field ({xj^i). The truncated expectation value is zero unless both the total number of 
fluctuation fields is even, and each cluster contains at least one fluctuation field. It is 
evaluated with Wick's theorem as a sum over pairings. The normal ordering forbids 
pairings within clusters. Step two is a re-normal ordering of the background fields. 
This is done with the help of 



j=l iG/j 



X 



■v-r 



E E Illls^i 



X 



■v-r 



JlC/i 



j=l 



n ■ n <pi^u. 

\j = l i&J'r 



(90) 



The Gaussian expectation value is here un-truncated and with negative covariance. It 
is defined by Wick's rule as a sum of all pairings. The number of fields being contracted 
is again required to be even. The expectation value is zero for an odd number of fields. 
The sums in (^) go over all subsets Jj of Ij including the empty set and Ij itself. Their 
meaning is that for each field in cluster j we decide whether it be a truly external field 
or contracted with a field in another cluster j' ^ j. We can think of this process as a 
generation of additional loops with normal ordering covariance in every term generated 
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by the fluctuation integral. Having done both the truncated fluctuation integral and 
the un-truncated normal ordering "integral" we obtain 



/ r,<s?/' 



E 

p=i 



1 



yi 



•d 



D 



y2p 



{2p)\ 

''-\Ly,,...,Ly,p) 



(91) 



with kernels given by 



^2p 



(y„...,yg = (2p)! E ■•• E E ■■■ E 

/lC{l,...,2fcl} I,nC{l,...,2k,n} JlQh JmCI„ 



Vl| + -+I^m|,2p 



/ f ft n ft (7^ 

\j=ii€{i,...,2kj}\.j, y j=i Vl^^-jJ- 



ie{i,...,2kj}\ij 



T 



(92) 



The external points on the right side are here understood to be renamed. There order 
is of no importance since all kernels are symmetric. The terms on the right hand side 
depend on the cardinalities and \Jj\ only. We are now ready to do the estimates. 
The right hand side of (021) is a sum of terms coming from writing out explicitely the 
expectation values in terms of contractions, the Feynman amplitudes. Their precise 
form is of no importance here. We restrict our attention to the number of lines and 
loops and combinatoric factors and find 



(93) 



L = 


2E(2^- 




Lh = 


-| m 

^E(2A:,- 




T tree 

Lifj — 


m — 1, 




Ls = 


1 m 

^Ed^.l- 





(94) 

(95) 
(96) 



Here L is the total number of lines, that is, factors of propagators. Lh is the number 
of hard lines of fluctuation propagators coming from the truncated expectation value. 
Ls is the number of soft lines from re-normal ordering. Each contraction is cluster 
connected in terms of hard lines. Therefore each contraction contains a tree of hard 
lines. This number of tree hard lines is L^*^^. The difference of L and L^*^^ is the 
number of loop integrals. We are ready to do the estimate. Put Li _2e-norms on all 
loop lines of hard and soft propagators. Put Loo,-2e-norms on the hard tree lines. Put 
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Loo,e-norms on all vertices. Each of the terms on the right side of ( p2D being bounded 
we find 

_ 2fc,-l 2fe„-l |/l| 

\\Kl;-'-^\U,<{2p)\ E ■■■ E E ••• E hA+-MM,2r> 

|/l|=0 |7„| = l|Ji|=0 |Jm|=0 

n /^/^ W \Lj\ \ 1 r^a,) \ iipnir^ ||fll^«"''^''l7;ll^s 

M vv i^^-i / V i-^^-i / 2^;)! " "°°'7 " "''"'^ " "''"'^ 

(2Lh-1)!!(2L5-1)!!. (97) 

More sophisticated estimates will be presented elsewhere. Hence ||-f^2p ' '''"''lloo,e is finite 
in momentum space. Furthermore all loop integrals converge, and the result is smooth 
in the external momenta. It can in particular be Taylor expanded to any desired order. 
But then we are done: the right hand side of the recursion relation is smooth and can 
therefore be solved as above. Thus the expansion is indeed finite to every order of 
perturbation theory. 



9 Second order 

The second order scaling equation ( |39D involves a truncated expectation value of two 
first order interactions. The first order interaction is given by eq. (0), including a 
0^-vertex and a wave function term. The first order wave function constant will come 
out of this second order computation. The expectation value is 

{V^'^; V^'^f^^^^ = (O4,o; O,,o)lsi. + 2Ci {O,,o; 02,2)lsi. + Ci (^2,2; 02,2)ls^ . (98) 

Each contribution is best taken care of separately. The computation consists of three 
steps. The first step is the fluctuation integral. It generates contractions between the 
vertices and shifts the normal ordering covariance. The second step is the rescaling 
of the external field. It restores the invariant normal ordering covariance and rescales 
the kernels of the effective interactions. The third step is a re-normal ordering of the 
result. Re-normal ordering creates additional loops with normal ordering propagators. 
The final result is 

(04,0; o,,o)ls^ = IJ^^'^J ^''y ■■ i^i^My) ■ l'^"" (^r(Lx - Lyf + 

T{Lx — LyYu{Lx — Ly) + V{Lx — Ly)u{Lx — LyY^ + 
l/d'^xj d^y : HxfijiyY : L'[3T{Lx - Lyf + 
6r(Lx — Ly)u{Lx — Ly) ) + 
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i |d^x I d^y : ^l^ixf^l^iyf : 20L6-^r(Lx - Ly), (99) 



2{0,,o;02,2fr,s^ = \jd^x:^{xf:L'i^jd^yT{y)^T{y) + 



+ 



2 / d^y.(i/)^r(y)^ 

I fd^'xl d^y : : 9>L^^T{Lx - Ly), (100) 

Here u{x) = L~'^^^v{L~^x) denotes a rescaled normal ordering covariance. The par- 
tial derivatives are meant to act on rescaled fluctuation propagators. The relevant 
part can be extracted by Taylor expansion of the Fourier transformed kernels. The 
inhomogeneous terms of the second order scaling equation are given by 



(pi,P2) = L-^^^-^) [hCiPl + L\M + 2c,) + L^-^Ci'f (L-Vi) + 

—V * f ★ f(L- Vi) + ^^f ★ f * u{L-^pi) + L'^ricuic u{L-^pM , (102) 

(Pi, ■ ■ ■ ,P4) = L-'('-^^ (b2 + SL'-^plr{L-'p,) + 

3L^-^f ^ f(L- Vi + L-^P2) + GL^'^f ^ u{L-^pi + L-^p,)^ (103) 
(Pi, ■ ■ ■ ,P6) = L-'^'-''^ (20L«-^^f(L-Vi + L-^p2 + L-'p,)) , (104) 

symmetrized in their momenta. The sum of momenta is constrained to zero. Here -k 
denotes convolution in momentum space (times {2n)~^). The constants stand for the 
convergent loop integrals 



= /7T^4p'^(p)^ (105) 



2p/ N2 

(2^^ 

C2 = / 7^/r(p)S(p). (106) 



(2 



TT 



The momentum space convolutions are conveniently computed with the help of the 
parameter representations 

r(p) = / dae-'^P (107) 

m(p) = / dae-^'P . (108) 
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The momentum space integrals become Gaussian and can be reduced to the two for- 
mulas 

/ (^''~''''''~"'^"'^' " (47r)^(Q;i + a2)^ei^^', (109) 

gcija2 + aia3 + c<2a3 ^ _ (110) 

A welcome feature of this parameter representation is that the result depends only 

on the external momentum squared. The scaling equation for the coefficient 62 in the 
/3-f unction follows from evaluation at zero momentum. It is given by 

62 = (3f*f(0) + 6f*S(0)) . (Ill) 

Its four dimensional value with exponential cutoff is computed to 

, _ -61og(L) 



{Any 



(112) 



It follows that the second order flow on the renormalized trajectory in four dimensions 
is given by 

P{9)-9-'-^9' + 0(g^). (113) 



(47r)2 

It follows that the model is asymptotically free in the infrared limit at weak coupling. 
The scaling equation for the first order wave function constant is found to be 



1 d f 1 

Cl = -^-qTt; *f *f(p)+f*f*S(p)+f*{2^M(p) 



62 %') Vs 

The result in four dimensions is 



(114) 

p2=0 



There is one relevant coordinate left at second order, the quadratic kernel at zero 
momentum. It is directly determined by the scaling equation 



(p--^^^-^) -^V?\^^^) ^ L-2(^-^)(l2Ci(ci + 2c2) + 



—V ★ f ★ f(0) + L^f ★ f ★ S(0) + L^f ★ S ★ S(0Ol6) 
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Doing these integrals as a last exercise the effective mass constant at second order in 
four dimensions turns out as 

^^'^(0,0) = ^ (21og(2) -log(3) - 1) . (117) 

This completes the computation of the second order renormalization constants. The full 
second order kernels are then given by summed scaled parameter integrals. The explicit 
computation of this irrelevant part will not be pursued further here. If this second 
order approximation is for instance used in a numerical simulation as renormalization 
group improved action, it is reasonable to approximate it by a momentum space Taylor 
expansion to some chosen degree of irrelevancy. Then this computation supplies the 
non- irrelevant part. 



10 Conclusions 

The standard renormalization scheme departs from a bare action. The renormalized 
trajectory is reached upon infinite iteration of block spin transformations. A further 
block spin transformation does no harm to a renormalized action. It merely generates 
a renormalized renormalization group flow on the renormalized trajectory. In view of 
this renormalization group flow on the renormalized trajectory we speak of running 
couplings. We should mention that the behavior of a renormalized action under the 
field theoretic renormalization group was used by Callan [ |C76| | in his proof of the BPHZ 



theorem, a polished version of which was given by Lesniewski ||L83|] . By universality the 
trajectory can be approached from a variety of bare actions. The notion of bare action 
is not unique. Consider for instance again the ultraviolet limit of (/)^-theory in four 
dimensions (say at negative coupling). The bare action is required to converge to the 
stable manifold of the trivial fixed point. But this stable manifold of irrelevant couplings 
is infinite dimensional. The approach presented here is free of bare ambiguities. It is 
designed upon equations for the renormalized action. 

The renormalized trajectory can be viewed itself as a means to investigate ultravi- 
olet and infrared limit of a Euclidean field theory. It is ideally suited to perform an 
infinite number of block spin transformations. A single block spin step on the trajectory 
translates to a (generally non-linear) transformation on the low dimensional space of 
running couplings. As a finite dimensional dynamical system it is comparatively easy 
to analyze. Consider for instance the 0'^-trajectory in the formalism above. There a 
block spin transformation becomes a transformation of a running coupling g in terms of 
a step /9-function (3{g) = big+b2g^ + 0{g^) with coefficients bi = L^~^ and 62 < 0. Here 
L > 1 is the block scale and D is the Euclidean dimension. For D = 4 it follows that 
the flow on the 0'^-trajectory is not asymptotically free at weak positive coupling in the 
ultraviolet limit. The running coupling shrinks under a block spin transformation. An 
infinite number of block spin transformations require either the renormalized coupling 
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to approach zero or the bare couphng (as is beheved) to tend to infinity. Perturba- 
tion theory is in this sense renormahzation group inconsistent for the ultraviolet limit 
in four dimensions. We remark that another non-trivial renormahzation group fixed 
point on the renormalized trajectory would allow to perform the ultraviolet limit. It 
would albeit suffice to require an infinite number of inverse block spin transformations 
for the running coupling to diverge. Either way requires control of the renormalized 
trajectory outside a vicinity of the trivial fixed point. This question is believed to be 
settled in favour of triviality by Aizenmann and Frohlich. See |[FFS92|| and references 
therein. However even in this situation of an inconsistency of perturbation theory in 
the ultraviolet limit perturbation theory for the 0^-trajectory at weak couplings is a 
well posed problem. 

The computational scheme presented here is very fiexible. It applies to any form 
of block spin transformations in the vicinity of a trivial fixed point. The expansion is 
iterative and built upon a fairly simple recursion relation. One can therefore expect 
it to be useful for accurate bounds on renormalized series. For instance tree decay of 
interactions is an immediate consequence. However the Loo-bound is already sharper 
than the estimates of Polchinski [ [F84|| , and the improvements by Keller, Kopper, and 
Salmhofer ||KKS90|| , and by Hurd |P389|| . It does not lead to logarithms in a renormal- 
ization scale. Other sophisticated expansion technology can be found in Rivasseau's 
textbook [[R91|| together with a set of references to the Paris school of renormahzation. 
In particular Rivasseau's effective expansion seems to be related to the one considered 
here. 

The question whether our approach makes sense beyond perturbation theory is 
a very interesting one. Whether renormahzation group invariance plus certain intial 
conditions determine a renormalized trajectory beyond perturbation theory we do not 
know. It is certainly the case where the perturbation expansion converges. Last not 
least it would be highly desirable to develop iterative techniques for unstable manifolds 
of nontrivial fixed points. A clear presentation of the weak coupling case around the 
trivial fixed point could be part of this way. 
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